Perturbative interpretation of relativistic symmetries in nuclei 
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Perturbation theory is used systematically to investigate the symmetries of the Dirac Hamilto- 
nian and their breaking in atomic nuclei. Using the perturbation corrections to the single-particle 
energies and wave functions, the link between the single-particle states in realistic nuclei and their 
counterparts in the symmetry limits is discussed. It is shown that the limit of S — V = const and 
relativistic harmonic oscillator (RHO) potentials can be connected to the actual Dirac Hamiltonian 
by the perturbation method, while the limit of S + V = const cannot, where S and V are the 
scalar and vector potentials, respectively. This indicates that the realistic system can be treated as 
a perturbation of spin-symmetric Hamiltonians, and the energy splitting of the pseudospin doublets 
can be regarded as a result of small perturbation around the Hamiltonian with RHO potentials, 
where the pseudospin doublets are quasidegenerate. 
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It is well known that the spin symmetry (SS) breaking, 
i.e., the remarkable spin-orbit splitting for the spin dou- 
blets (n, l,j = I ± 1 /2), is one of the most important con- 
cepts for understanding the traditional magic numbers 
(2, 8, 20, 28, ...) in atomic nuclei 0, Meanwhile, a 
new symmetry, the so-called pseudospin symmetry (PSS) 
0, 13] , is introduced to explain the near degeneracy be- 
tween two single-particle states with the quantum num- 
bers (n-l,l + 2,j = I + 3/2) and (n,l,j = I + 1/2) by 
defining the pseudospin doublets (n — n—l,l = l + = 
I ±1/2). The splittings of both spin and pseudospin dou- 
blets play critical roles in the shell structure evolutions. 
Thus, it is a fundamental task to explore the origin of SS 
and PSS, as well as the mechanism of their breaking. 

Since the suggestion of PSS in atomic nuclei, there 
have been comprehensive efforts to understand its ori- 
gin. Apart from the rather formal relabeling of quantum 
numbers, various proposals for an explicit transformation 
from the normal scheme to the pseudospin scheme have 
been discussed [B-Q • Based on the single- particle Hamil- 
tonian of the oscillator shell model, the origin of PSS is 
connected with the special ratio in the strength of the 
spin-orbit and orbit-orbit interactions The relation 
between the PSS and the relativistic mean field (RMF) 
theory Q was first noted in Ref. in which Bahri et 
al. found that the RMF theory approximately explains 
such a special ratio in the strength of the spin-orbit and 
orbit-orbit interactions. 

As substantial progress, the PSS was shown to be a 
symmetry of the Dirac Hamiltonian, where the pseudo- 
orbital angular momentum I is nothing but the orbital 
angular momentum of the lower component of the Dirac 
spinor, and the equality in magnitude but difference in 
sign of the scalar potential 5(r) and vector potential V(r) 
was suggested as the exact PSS limit by reducing the 
Dirac equation to a Schrodinger-like equation fioj . As a 



more general condition, d(S + V)/dr = can be approxi- 
mately satisfied in exotic nuclei with highly diffuse poten- 
tials jll|, [ljj]- Meanwhile, based on this limit, the pseu- 
dospin SU(2) algebra was established [l3j . and the spe- 
cific node structures of the pseudospin doublets were illu- 
minated (l4| . Furthermore, the Dirac Hamiltonian with 
spin and pseudospin SU(2) symmetries can also be de- 
rived in supersymmetric (SUSY) patterns [l5j . However, 
since there exist no bound nuclei within S + V = const, 
the non-perturbative nature of PSS in realistic nuclei has 
been presented in Refs. which is also related to 

the consideration of the PSS as being a dynamical sym- 
metry [19( 1 . 

On the other hand, the relativistic harmonic oscilla- 
tor (RHO) potentials were used to understand the origin 
of PSS Ho, [ll|. Subsequently, the spin and pseudospin 
U(3) algebra was established in the Dirac Hamiltonian 
with RHO potentials [HI, [22j|. Recently, Typel pointed 
out that the Hamiltonian with spin U(3) symmetry is one 
of the simplest cases where the pseudospin symmetry- 
breaking potential derived in the SUSY framework van- 
ishes [23| . Meanwhile, Marcos et al. commented that the 
quasi-degeneracy of the pseudospin doublets in realistic 
nuclei can be considered as the breaking of their degener- 
acy in the Dirac Hamiltonian with RHO potentials 24 1 . 



In this Rapid Communication, the perturbation theory 
will be used for the first time to investigate the symme- 
tries of the Dirac Hamiltonian and their breaking in re- 
alistic nuclei. The perturbation corrections to the single- 
particle energies and wave functions will be accurately 
calculated numerically. In this way, the link between the 
single-particle states in realistic nuclei and their counter- 
parts in the symmetry limits will be constructed explic- 
itly. 

Assuming the spherical symmetry, the radial Dirac 
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FIG. 1: (Color online) Single-particle potentials for neutrons 
in the nucleus 132 Sn. The self-consistent potentials calculated 
by the RMF theory with PK1 |25( | are shown as solid lines. 
The potentials -A - M in Hq S and E + M in Hq ss are 
illustrated as dashed, and dash-dotted lines, respectively. 



equations can be cast in the form of 
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where E(r) = S(r) + V(r) and A(r) = S(r) - V(r) de- 
note the combinations of the scalar and vector potentials, 
and k is defined as k = (I — j)(2j + 1). Taking the nu- 
cleus 132 Sn as an example, the potentials E(r) and A(r) 
for neutrons calculated by the self-consistent RMF the- 
ory with the effective interaction PK1 (2f| are shown in 
Fig. [U It is generally found that E(r) ~ 70 MeV and 
A(r) ~ 700 MeV in the realistic nuclei. 

Following the idea of Rayleigh-Schrodinger perturba- 
tion theory, the Dirac Hamiltonian H in Eq. @ can be 
split as 



H = Hn 



w, 



or equivalcntly 



H = H - W, 



(3) 



(4) 



where Hq leads to the exact spin (pseudospin) symme- 
try and W is identified as the corresponding symmetry- 
breaking potential. The condition 
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where W m k = (W m ) W (^fc), determines whether W can 
be treated as a small perturbation and governs the con- 
vergence of the perturbation series prjj . 



In the case of the spin and pseudospin SU(2) symmetry 
limits shown in Ref. [l5j . the Dirac Hamiltonian with 
exact symmetries reads 
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whose eigenenergies are denoted as Eq in general, and 
the corresponding symmetry-breaking potentials are 
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In contrast to adopting the Schrodinger-like equations 
in the previous studies [l(| E3, [H, H3, it is clearly 
shown that the operators H, Hq and W used in the 
present calculations are all Hcrmitian, and they do not 
contain any singularity. This allows us to perform the 
order-by-order perturbation calculations. In addition, it 
should also be noticed that only W corresponds to the 
symmetry-breaking potential within the present decom- 
position, thus the ambiguity caused by the strong cancel- 
lations among the different terms in the Schrodinger-like 
equations can also be avoided. Therefore, the present 
method can provide a clear and quantitative way for in- 
vestigating the perturbative nature of SS and PSS. This 
method can be universally applied to the cases that the 
nature of the symmetry is either perturbative or non- 
perturbative. For the case where the nature of the sym- 
metry is perturbative, the link between the single-particle 
states in realistic nuclei and their counterparts in the 
symmetry limits can be constructed. For the case where 
the nature of the symmetry is non-perturbativc, the di- 
vergence of the perturbation series can be found explic- 
itly. 

In the present calculations, as illustrated with dashed 
and dash-dotted lines in Fig.[TJ the constant potentials in 
Eqs. © are chosen as -A - M = -350 MeV and S + 
M = 900 MeV. We have checked that the convergence 
of the perturbation calculations is not sensitive to these 
values. 

In Fig. O taking the spin doublets k = If and the 
pseudospin doublets k = Id as examples, the values of 
\W m k/(E m — Ek)\ are plotted as functions of the energy 
differences E m — Ek- For the SS case, the unperturbed 
cigenstates are chosen as those of Hq S in panel (a), while 
the unperturbed cigenstates are chosen as those of H in 
panel (b), namely, the former perturbation calculations 
are performed from a spin-symmetric Hamiltonian Hq S 
to the realistic Hamiltonian H , whereas the latter ones 
are performed from H to Hy S . For the PSS case, since 
there are no bound states in the pseudospin-symmetric 
Hamiltonian Hg SS , the perturbation calculations are per- 
formed only from H to Hg SS as shown in panel (c). For 
the completeness of the single-particle basis, the single- 
particle states m must include not only the states in the 
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FIG. 2: (Color online) Values of |FF mfc /(.E m - E u )\ vs the 
energy differences E m — E k for the spin doublets k — If 
(panels (a) and (b)) and the pseudospin doublets k — Id 
(panel (c)). The unperturbed eigenstates are chosen as those 
of Hq S (panel (a)) and H (panels (b) and (c)), respectively. 
The single-particle states m include the states in the Dirac 
sea and Fermi sea. 



Fermi sea, but also those in the Dirac sea. Since the 
spherical symmetry is adopted, only the states m and k 
with the same quantum numbers I and j lead to non- 
vanishing matrix elements W m k- It is seen that the val- 
ues of \W m k/(E m — Ef.)\ decrease as a general tendency 
when the energy differences \E m — increase. From 
the mathematical point of view, this property provides 
natural cut-offs of the single-particle states in the per- 
turbation calculations. 

Although the potentials satisfy |A — A| » |S — S |, 
the largest value of \W m k/(E m — Ek)\ is roughly 0.10 
(H$ s to H) or 0.06 (H to H$ s ) for the SS case, whereas 
it is about 0.6 for the PSS case because different compo- 
nents of the Dirac spinors are involved: 
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FIG. 3: (Color online) Single-particle energies of spin doublets 
If (panels (a) and (b)) and pseudospin doublets Id (panel 
(c)) obtained at the exact symmetry limits, and by the first-, 
second-, and third-order perturbation calculations, as well as 
those by the RMF theory. The unperturbed eigenstates are 
chosen as those of Hq S (panel (a)) and H (panels (b) and 
(c)), respectively. 



where for the Fermi states the upper component G(r) ~ 
0(1), and the lower component F(r) ~ 0(1/10). This 
indicates that the criterion in Eq. can be well fulfilled 
for the SS case, but questionable for the PSS case. 

Let us then examine the perturbation corrections to 
the single-particle energies of the spin doublets If and 
pseudospin doublets Id. In panel (a) of Fig. [3j by choos- 
ing the unperturbed eigenstates as those of Hq S , the 
single-particle energies obtained at the exact spin sym- 
metry limit, and their counterparts obtained by the first- 
, second-, and third-order perturbation calculations, as 
well as those obtained by the self-consistent RMF the- 
ory, are shown from left to right. Meanwhile, the corre- 
sponding results obtained by choosing the unperturbed 
eigenstates as those of H are shown in panels (b) and (c) 
of Fig. [3] It can be seen clearly that the spin-orbit split- 
ting is well reproduced by the second-order perturbation 
calculations as shown in panel (a), and in the reversed 
way the energy degeneracy of the spin doublets can be 
well restored as shown in panel (b). This verifies that 
for studying the relationship between the eigenstates of 
Hq and H by perturbation theory, it is equivalent to use 
the definitions H = H + W and H = H - W. In 
contrast, as shown in panel (c), the energy degeneracy 
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FIG. 4: (Color online) Upper panels: Same as Fig. [2] but for the case of the RHO potentials. Lower panels: Same as Fig. [3] 
but for all single-particle states in the pf major shell. The unperturbed eigenstates are chosen as those of Hq* 10 (panels (a) 
and (c)) and H (panels (b) and (d)), respectively. 



of the pseudospin doublets cannot be restored up to the 
third-order perturbation calculations, and there exist no 
bound eigenstates of Hq°. Thus, the link between the 
pseudospin doublets in realistic nuclei and their counter- 
parts in the S + V = const limit is unclear. In addi- 
tion, choosing the eigenstates of H as the unperturbed 
eigenstates, the perturbation corrections to the single- 
particle wave functions arc also evaluated. It is found 
that the identity of the wave functions for the spin dou- 
blets, Go(lf 5 / 2 ) = Go(l/7/2), can be well reproduced by 
the second-order perturbation corrections, but the iden- 
tity of the wave functions for the pseudospin doublets 
^0(1/5/2) = ^o(2p 3 / 2 ), cannot be fulfilled. 

Therefore, from the perturbative point of view, the 
bridge can be constructed to connect the Dirac Hamilto- 
nian in realistic nuclei with the symmetry limit of S— V = 
const, but not S + V = const. This indicates that the 
realistic system can be treated as a perturbation of the 
spin-symmetric Hamiltonian. This also confirms in an 
explicit way that the nature of PSS is non-perturbative, 
as the Dirac Hamiltonian with S + V = const is regarded 
as the symmetry limit. 

However, it has been pointed out that the energy split- 
ting of the pseudospin doublets in realistic nuclei could 
be alternatively considered as the breaking of their de- 
generacy appearing in the spin-symmetric Hamiltonian 
with RHO potentials p3l . [24| . In the following, we assess 
this statement in a perturbative way. 

In the case of the spin U(3) symmetry limit shown in 
Ref. [22|, the Dirac Hamiltonian H in Eq. ^ is split as 

H = H™° + W RHO , (9) 



with 

° "I f + 7 -Ao-Mj' (10) 

and 

W hho = (S-Eho Aq 0_ a ) 5 (n) 

where EhoM = c o + c 2 r2 has the form of a harmonic 
oscillator. Here, H^ HO leads to the energy degeneracy of 
the whole major shell, and IV RHO is identified as the cor- 
responding symmetry-breaking potential. In the present 
investigation, we choose —Aq—AI = —350 MeV as in Hq S 
and Co + M = 865 MeV. As discussed before, the pertur- 
bative properties are not sensitive to these two constants. 
Meanwhile, the coefficient C2 is chosen as 1.00 MeV/fm 2 
to minimize the perturbations to the pf states. 

In the upper panels of Fig. 21 the values of 
\W^°/(E m — Ek)\ for the pseudospin doublets k = Id 
are shown as functions of the energy differences E m — Ek ■ 
It is found that the general patterns shown in panels 

(a) and (b) are the same as those in panels (a) and 

(b) of Fig. [21 respectively, and the largest perturbation 
correction is roughly 0.16 (H^ HO to H) or 0.10 (H to 
Hq HO ). This indicates that the criterion in Eq. ([5]) is 
fulfilled, even though not as well as in the SS case. In the 
lower panels of Fig. 2J the single-particle energies of the 
states in the pf major shell obtained at the exact sym- 
metry limit and by the self-consistent RMF theory, as 
well as their counterparts obtained in the first-, second-, 
and third-order perturbation calculations, are shown. As 
shown in panel (c), not only the spin-orbit splitting, but 
also the pscudospin-orbit splitting, is well reproduced by 
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the third-order perturbation calculations, and in the re- 
versed way the energy degeneracy of all the states in the 
pf major shell can be well restored as shown in panel 
(d). Thus, the link between the pf states in realistic nu- 
clei and their counterparts in the symmetry limit with 
RHO potential can be explicitly established. Further- 
more, it is found that the single-particle wave functions 
of (ir RHO ) can also be reproduced by the second-order 
perturbation calculations starting from if RHO (H). 

Therefore, the bridge connecting the Dirac Hamilto- 
nian in realistic nuclei and that with RHO potentials 
can be clearly constructed using perturbation theory. 
This indicates that the energy splitting of the pseudospin 
doublets can be regarded as a result of small perturba- 
tion around the Dirac Hamiltonian with RHO poten- 
tials, where the degeneracy of the pseudospin doublets 
appears. 

Of course, one must always keep in mind that the ac- 
tual picture in nuclear spectra is generally more complex 
than that of a simple potential model. Still, there are 
well-identified cases where the core polarization effects 
are small enough (correlation and polarization diagrams 
may compensate each other and give relatively large spec- 
troscopic factors) to allow for the notion of single-particle 
state to hold. 

In summary, the symmetries of the Dirac Hamilto- 
nian and their breaking in realistic nuclei are investigated 
in the framework of perturbation theory. The present 
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